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Abstract 

A Dalitz-plot analysis of B^{t) —^pn-^ tt+tt^tt^ decays allows one to obtain the 
CP- violating phase a. In addition, one can extract the various tree (T) and penguin 
(P) amplitudes contributing to these decays. By comparing the measured value of 
\P/T\ with the theoretical prediction, one can detect the presence of physics beyond 
the standard model. 
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A great many methods have been proposed for obtaining information about the 
CP phases a, (3 and 7 of the unitarity triangle Almost all of these involve CP- 
violating asymmetries in hadronic B decays [2]. The aim is to test the standard 
model (SM) explanation of CP violation, and hopefully find evidence for physics 
beyond the SM. 

The cleanest methods (i.e. those in which the theoretical hadronic uncertainties 
are very small) involve B decays which are dominated by a single amplitude, such 
as -B^(t) J/ipKg. However, many decays receive contributions from both tree and 
penguin diagrams with different weak phases [S]. A-priori, one would think that one 
cannot obtain clean phase information from such decays. Fortunately, techniques 
have been developed for removing the unwanted "penguin pollution." For example, 
an isospin analysis oi B tttt decays allows one to remove this contamination and 
obtain sin2Q; cleanly f4j, albeit with discrete ambiguities. 

In fact, this isospin analysis gives us even more information. In particular, one 
can also obtain the magnitudes and relative phases of the tree (T) and penguin (P) 
amplitudes in B^ vr+vr" It is therefore possible to compare the experimental 
value of \P/T\ with that predicted by theory. If a significant discrepancy is observed, 
it would signal new physics p^. 

An alternative technique for obtaining a involves B ^ pn decays. By performing 
a Dalitz-plot analysis of -B°(t) —>■ pn —>■ 7r+7r~7r° decays, one can remove the penguin 
contributions from B —>■ pn decays and obtain a [7j. Compared to i? ^ vrvr, the 
advantage of this method is that it is possible to extract both sin2Q; and cos 2a, 
so that one obtains 2a with no discrete ambiguity. Another advantage is that 
it is not necessary to measure processes involving two final-state 7r° mesons. The 
disadvantage of this method is that one must understand the continuum background 
to such decays with considerable accuracy, as well as the correct description of 
p ^ 7T7T decays, and these may be difficult. 

Here too there is enough information to obtain the magnitudes and relative 
phases of the tree and penguin amplitudes. Thus, one can measure \P/T\ in B ^ pn. 
As in i? — s> TTvr, a comparison of this ratio with the theoretical prediction can reveal 
the presence of new physics. In this paper we perform such an analysis. As we will 
show, the B ^ pn method has two advantages compared to -B ^ vrvr for searching 
for physics beyond the SM in this way. First, the fact that there is no discrete 
ambiguity in 2a improves the prospects for finding new physics. Second, the \P/T\ 
ratio is expected to be smaller than in B tttt, which makes it easier to see a 
new-physics signal, should it be present. 

We begin with a brief review of the B ^ pn Dalitz-plot analysis within the SM 
[7]. There are five B —>■ pn amplitudes which satisfy a pentagon isospin relation. 
All amplitudes receive contributions from both tree and b —>■ d penguin amplitudes. 
The tree amplitude is proportional to V*i,Vud, while the penguin amplitude has con- 
tributions from internal u, c and t quarks, proportional to V*ijVud, Vetoed and V^^V^d, 
respectively. Using the unitarity of the CKM matrix, Vf*^Vtb + V*^Vcb + KdKb = 0, 
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we can eliminate the c-quark contribution. Furthermore, the piece proportional to 
Vub^ud can be absorbed into the tree amplitude. Thus, the penguin amplitude in- 
cludes only the piece proportional to V^*^Vtb. It is convenient to rescale the amplitudes 
by e*^, leading to the following expressions for the amplitudes: 



'S'+o 






) = 


y+Og-ia _^ p+0 ^ 


So+ 


= e"^V2A{B^ 


- ^ p%+: 


) = 


rpO+^-ia ^ pO+ ^ 


S+- 


= e*^A(5°^ 












p-7l + ) = 


T- 


-+e-'" + P-+ , 


Soo 




. p\') = 




-0 + 'S'o+ — — — 



where = — P+°. In the above we have explicitly written the weak phase a, 
while the Tj and the Pi include strong phases. (Throughout the paper, we use the 
subscript 'i' to denote all of the pvr charge combinations: i = +0, 0+, H — , — h, 00.) 
The corresponding amplitudes for the CP-conjugate processes. Si, are obtained by 
changing the signs of the weak phases. 

The key point is that all of the neutral — > pvr amplitudes contribute to 
B^ 7r"'"7r~7r°. We can therefore write 

A(P0 ^ tt+tt-ttO) = + + fSoo/2 , (2) 

where the are the kinematic distribution functions for the pions produced in the 
decay of the p*. The B^ mesons can decay to the same final state: 

A(P0 tt+tt-ttO) = f-S+. + f+S.+ + fSoo/2 . (3) 

The time-dependent measurement of the Dalitz plot for B^{t) — > vr'^Tr^'vr^ then allows 
one to extract the magnitudes and relative phases of each of the Sij and Sij in 
Eqs. (121) and © IZj. By taking the ratio of the relations 

S+. + + Soo = iT+'^ + T°+) e-*" , 

S+- + S-+ + S00 = (T+o + T0+) e'" , (4) 

one obtains e~^*". We therefore see that, using this method, the CP phase 2a can 
be extracted with no ambiguity. 

It is also possible to obtain the tree and penguin contributions to the amplitudes 
in Eq. ((T}. We define the following observables: 
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Here Bi, ai and "io-lff are, respectively, the branching ratio, direct CP asymmetry, 
and measure of indirect CP violation for each decay. We remark that each of 2a^jj 
and 2agjj are automatically zero since they involve charged B decays. (Note: the 
indirect CP asymmetry is usually written with an explicit mixing phase q/p = e~^*^. 
This phase is removed when one rescales the amplitudes by e*^ as in Eq. (Q).) We 
have 



— = —2i sin a Ti 
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Si.e''' - S.e-'" = -2i sin a Pi. (6) 
It is then straightfoward to obtain |Tjp and l-Pjp: 



\Ti\ — Ri 
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where 



Ri = 
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Note that Ri is proportional to B^ [Eq. (0)]. The proportionality constant depends 
on which decay is being considered, see Eq. ((H). 

Suppose now that there is physics beyond the SM. If present, it will affect mainly 
B^-B^ mixing and/or the b ^ d penguin amplitude. In the SM, the weak phase of 
B^-B^ mixing is equal to that of the t-quark contribution to the b ^ d penguin. 
This is reflected in the fact that the weak phase multiplying the term Pj in Eq. (Q) 
is zero. If new physics is present, these two weak phases may be different. One 
can take this possibility into account by including a new-physics phase O^p in the 
B ^ pn amplitudes: 
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(9) 

The extraction of a is unchanged by the presence of the new-physics parameter ^jvp 
(though its value may include new contributions to P^-P^ mixing). However, the 
expressions for Tj and Pi are modified. We now have 

Si e^^^^ - Si e"^^^^ = -2i sin (a - O^p) Ti , 

Si e^" - Si e"*" = -2i sin(a - O^p) Pi , (10) 
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so that 
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1 - cos(2a - 20jvp) 
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(11) 



The expressions for the Tj and Pj are therefore ahered in the presence of new physics. 
Thus, by comparing the measured value of a particular \P/T\ with that predicted 
by theory (within the SM), we can detect the presence of a nonzero 6np- (It is also 
possible for new physics to affect the magnitudes of the Tj and Pj. This possibility 
is implicitly included in our method.) 

The first step is therefore to evaluate the theoretical value of \P/T\. However, 
there are many \P/T\ ratios that can be considered. We concentrate only on the 
(color-allowed) neutral decays P^ p^vr^. There are several reasons for this. First, 
the Dalitz plots for the charged B decays are much more difficult to obtain since 
they require the detection of two vr^'s. Second, the branching ratio for the color- 
suppressed decay — > p°7r° is expected to be quite a bit smaller than those of 



p TT^. Finally, below we will use QCD factorization to estimate the theoretical 
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size of the |P/T| ratios, and nonfactorizable effects are expected to be small for 
color-allowed decays. 

The value of |P/T| for P^ p^vr^ has been calculated in the literature. In 
Ref. [H], this was done using naive factorization and including only the t-quark 
contribution to the b ^ d penguin. A more recent computation has been done 
by Beneke and Neubert [H] in the context of QCD factorization FTUl. Since QCD 
factorization is a state-of-the-art framework, using expansions in l/rrib and a^, we 
will follow this approach. (It should be noted, however, that the |P/T| ranges given 
in Refs. |H] and P are quite similar.) 
We define 
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Ref. P gives 
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and rus, and the 



The errors come principally from three sources: the values of \Vub 
size of weak annihilation effects. Note that the two ratios are determined by very 
different dynamics, so that their near equality is a numerical accident. 

It is now necessary to decide on the numerical ranges to use for r~^~ and r in 
the analysis. Since the goal is to search for physics beyond the SM, it is important 
to be as conservative as possible. With this in mind, we will take the theoretical 
ranges for r~^~ and r~~^ within the SM to be 



0.05 < < 0.25 , 0.05 < r~+ < 0.25 . 



(14) 
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The above ranges are larger than those given in Eq. particularly on the upper 
side. We note that QCD factorization cannot account for the observed i? — > vrvr 
and B —>■ pn branching ratios [9J. Assuming no new physics — and the analysis of 
this paper can be used to test for such effects — there must be some contribution 
which is larger than its QCD factorization value. The enlarged ranges of Eq. (fT^ 
take this into account, as well as potential underestimates of factorizable errors 
(e.g. electroweak-penguin effects) and nonfactorizable effects. With the ranges of 
Eq. (HH), a significant discrepancy between the measured value of \P/T\ and its SM 
prediction will clearly be a sign of new physics. That is, within the SM, we expect 



for z = H — , — h. If it is found that the observables do not respect this inequality, 
this points to the presence of physics beyond the SM. 

We note in passing that, while the \P/T\ range in — > pvr is ~ 10%, in B ^ tttt 
it is expected to be ~ 20-30% 13 HI] • We therefore conclude that the penguin 
pollution is likely to be more significant in 5 — > tttt Thus, if new physics is 
present, it will be easier to detect in _B — > pvr. 

As noted earlier, the CP phase 2a can be extracted from the B ^ pn method. 
However, this is not easy experimentally. In our analysis we therefore consider 
the possibility that only sin2Q; is measured (in which case one obtains 2a with a 
twofold ambiguity), as well as the case where 2a is known without ambiguity. In 
both scenarios, we consider two possible ranges for 2a: (i) 120° < 2a < 135°, (ii) 
165° < 2a < 180°, which can be considered to take into account the experimental 
errors in the measurements. 

Our results are shown in Fig. ^ We consider the two ranges for 2a given above. 
For each of these ranges. Fig. [T] shows the (correlated) allowed values of 2aljrj and 
Oj {i = H — , — h) that are consistent, within the SM, with the assumed range for 2a 
and the theoretical range for |PYT*| [Eq. (jl4j) ]. If only sin 2a has been measured, 
then, for a given range of 2a, both regions in Fig.^are allowed. If cos 2a can also be 
measured, then the left-hand region can be removed. In either of these scenarios, if 
the measured values of the observables do not lie within the SM region, this means 
that new physics — i.e. a nonzero 6,^p — is present. As can be seen from Fig. Q 
the new-physics region is quite large, so that we have a good chance of detecting 
the new physics via this method, should it be present. 

If no signal for new physics is detected, one can place an upper limit on the value 
of 6^p via 



To summarize, the measurement of the Dalitz plot of B^{t) ^ pn ^ tt+tt tt' 
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Figure 1: The region in 2a^f/-ai space {i = H — , — h) which is consistent with the 
theoretical prediction for \P^/T^\ [Eqs. fjl2|) . ()14|)]. for two ranges of 2a, shown above 
each figure. It is assumed that only sin 2a is measured, so that 2a is obtained with 
a twofold ambiguity. (This is the source of the two regions in both figures.) If cos 2a 
is also measured, allowing one to obtain 2a without ambiguity, the left-hand region 
must be removed in both figures. 



decays allows one to cleanly extract the CP-violating phase 2a, with no discrete 
ambiguity. One can also obtain the individual tree (T) and penguin (P) amplitudes 
in these decays. By comparing the measured value of a particular \P/T\ ratio with 
that predicted by theory, one can detect the presence of physics beyond the standard 
model. From both a theoretical and experimental point of view, the best \P/T\ ratios 
are those for p^vr^. The conservative ranges for these ratios are taken to be 

0.05 < |PYT*| < 0.25 (i = H — , — h). The region in (2agjj,aj) parameter space 
(2agjj and are, respectively, the measured indirect and direct CP asymmetries) 
which corresponds to this range of \P^/T^\ is relatively small. This therefore provides 
a good way of detecting the presence of new physics. 
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